Continuous-time quantum walks on one-dimension regular networks 
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In this paper, we consider continuous-time quantum walks (CTQWs) on one-dimension ring lattice 
of N nodes in which every node is connected to its 2m nearest neighbors (m on either side). In the 
framework of Bloch function ansatz, we calculate the spacetime transition probabilities between two 
nodes of the lattice. We find that the transport of CTQWs between two different nodes is faster 
than that of the classical continuous-time random walk (CTRWs). The transport speed, which is 
defined by the ratio of the shortest path length and propagating time, increases with the connectivity 
parameter m for both the CTQWs and CTRWs. For fixed parameter m, the transport of CTRWs 
gets slow with the increase of the shortest distance while the transport (speed) of CTQWs turns out 
to be a constant value. In the long time limit, depending on the network size N and connectivity 
parameter m, the limiting probability distributions of CTQWs show various patterns. When the 
network size N is an even number, the probability of being at the original node differs from that of 
being at the opposite node, which also depends on the precise value of parameter m. 

PACS numbers: 05.60.Gg, 03.67.-a, 05.40.-a 
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I. INTRODUCTION 



Quantum walks have important applications in various 
fields of solid-state physics, polymer chemistry, biology, 
astronomy, mathematics and computer science 0, [U, 0,0, 
A quantum random walk (QRW) is a natural exten- 
sion to the quantum world of the ubiquitous classical ran- 
dom walk. It was first introduced in [3] and extensively 
investigated recently in connection with possible applica- 
tions to quantum algorithms [8j]. The behavior of quan- 
tum walks differs from that of the classical random walks 
in several striking ways, due to the fact that quantum 
walks exhibit interference patterns whereas the classical 
random walks do not. For instance, the mixing times, 
hitting times and exit probabilities of quantum walks 
can differ significantly from analogously defined random 
walks 0, [To . [Tl| . In recent years, two types of quan- 
tum walks exist in the literature: the discrete-time quan- 
tum coined walks and continuous-time quantum walks 
[l2l [l3j . Although both the two types of quantum walks 
have efficient quantum algorithms with respect to their 
classical counterparts, quantum walks show some advan- 
tages in dealing with decoherence processes compared to 
the discrete-time quantum algorithms, which are very 
sensitive to environmental quantum noise 

Here, we focus on continuous-time quantum walks 
(CTQWs). Most of previous studies consider CTQWs 
on simple structures, such as, the line [H, cycle 
0,111, hypercube [lj], Cayley tree [I|, dendrimers [U 
and other regular networks with simple topology. Al- 
though CTQWs have received much attention and there 
has been some work about CTQWs on general graphs, 
many questions about CTQWs appear to be quite diffi- 
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cult to answer at the present time. For simple structures 
these quantum walks are analytical solvable and directly 
related to the well-known problems in solid state physics. 
Recently, Oliver Miilken et al have studied the space- 
time structures of CTQWs on one-dimensional and two- 
dimensional lattices with periodic boundary conditions 
[22l . The topology of the lattices they considered is 
oversimplified, i.e., each node is only connected to its two 
nearest neighbors. For regular graphs with symmetrical 
structure, the dynamics of the quantum transport is de- 
termined by the topology of the network. To this end, 
it is natural to consider quantum transport on general 
lattices with more connectivity. 



In this paper, we study CTQWs on one-dimension ring 
lattice of N nodes in which every node is connected to 
its 2m nearest neighbors (m on either side). This gener- 
alized regular network has broad applications in various 
coupled dynamical systems, including biological oscilla- 
tors [24j . Josephson junction arrays (2a |. neural networks 
[III, synchronization (27[, small- world networks (28| and 
many other self-organizing systems. We analyze quan- 
tum walks on such general network with periodic bound- 
ary conditions using the Bloch function approach [29j . 
which is commonly used in solid state physics. We de- 
rive analytical expressions for the transition probabilities 
between two nodes of the networks, and compare them 
with the results of CTRWs. 



The paper is structured as follows: In Sec. II we re- 
view the properties of CTQWs presented in Ref. [30( and 
give the exact solutions to the transition probabilities on 
the general ring network. Section III presents the time 
evolution of the probabilities. In Sec. IV, we consider the 
distributions of long time limiting probabilities. Conclu- 
sions and discussions are given in the last part, Sec. V. 
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II. CONTINUOUS-TIME QUANTUM WALKS 

Keeping in line with previous results on quantum 
walks, we study continuous-time quantum walk on net- 
works and compare the results with the classical coun- 
terparts. 



Continuous-time quantum walks on general 
networks 



We consider a walk on a general graph, which is a 
collection of connected nodes and simple links without 
weight and directions. The topology of such simple 
graphs can be described by the corresponding Laplace 
matrix A. The nondiagonal elements Aij equal to —1 if 
nodes i and j are connected and otherwise. The di- 
agonal elements An equal the number of total links con- 
nected to node i, i.e., An equals to the degree of node 
i. Classically, the evolution of continuous-time random 
walk is governed by the master equation Q 



dpk,j 
dt 



(1) 



Where Pk,j(t) is the conditional probability to find the 
CTRW at time t at node k when starting at node j. 
Matrix T is the transfer matrix of the walk, and relates 
to the Laplace matrix by T = — 7A Here, for the sake 
of simplicity, we assume the transmission rate 7 for all 
connections to be equal. Then the solution of the above 
equation is 



Pk,j{t) =< k\e u \j > 



(2) 



Quantum mechanically, the dynamical evolution equa- 
tion of continuous-time quantum walks is obtained by 
replacing the Hamiltonian of the system by the clas- 
sical transfer matrix, H = —T [{J [l(|. The states 
\j > endowed with the nodes j of the network form 
a complete, ortho-normalised basis set, which span the 
whole accessible Hilbert space, i.e., J2k \ k >< k \ = 1' 
< k\j >— Skj. The time evolution of state \j > is given 
by the Schrodinger Equation (SE) 



. d\j > 
1 dt 



H\j >, 



(3) 



Where the mass to = 1 and h = 1 is assumed in the above 
equation. Starting at time t from the state \j >, the evo- 
lution equation of the state \j > is \j,t >= U(t,t )\j >, 
where U(t,to) — e - l ^ f (*- t o) j s the quantum mechani- 
cal time evolution operator. The transition amplitude 
ctk,j(t) from state \ j > at time to state \k > at time t 



a kJ (t) =<k\e- itH \j>, 
Combining Eq. @, we have 
.da-kj 

1 



dt 



^Hkiaijit), 



(4) 



(5) 



We note that the different normalization for CTRWs and 
CTQWs. For CTRWs, Y, k Pk,j = 1 and quantum me- 
chanically J2k \ a k.j\ 2 = 1 holds. 

To get the exact solution of Eqs. |T|) and ([5]), all 
the eigenvalues and eigenvectors of the transfer opera- 
tor and Hamiltonian are required. We use E n to repre- 
sent the nth eigenvalue of A and denote the orthonor- 
malized eigenstate of Hamiltonian by \q n >, such that 
J2n \ln >< <7n I = I- The classical transition probability 
between two nodes is given by 



Pk,j 



(t)=Y. e ^ tEn < k \ln><q n \j >, 



(6) 



And the quantum mechanical transition probability be- 
tween k and j is 

7T*j(t) = \a Ko {t)\ 2 

= Y.nl^ lt(En ' El) < k\q n >< q n \j >< k\ qi X Ql \j > 

. . (?) 

For finite networks, TTkj(t) do not decay ad infinitum 
but at some time fluctuates about a constant value. This 
value is determined by the long time average of TTk,j(t) 

Xk,j = limr^oo i? Jq TTfc j (t)dt 

= En,l < k kn >< 1n\j >< k\qi >< qi\j > 

xlim T ^oc^/ T e- i ^ E «- E ')dt 
= E, t i ^E^.E, < k\q n >< q n \j >< k\qi >< qi\j > . 

(8) 



B. Continuous-time quantum walks on ID ring 
lattice and Bloch ansatz solutions 



In the subsequent calculation, we restrict our atten- 
tion on CTQWs on the general one-dimension ring lat- 
tices with periodic boundary conditions. The network 
organizes in a very regular manner, i.e., each node of the 
lattice is connected to its 2m nearest neighbors (to on ei- 
ther side), thus the Laplace matrix A takes the following 
form, 



2m, if i — j, 
An = i —1, if i = j ± z, z € [1, to] 
0, Otherwise. 



(9) 



The Hamiltonian of the system is given by H = 7 A For 
simplicity of analytical treatment, we set 7 = 1 in further 
calculations. The Hamiltonian acting on the state \j > 
can be written as 

m 

H\j >= (2m+ > - ^2 \j + z>, z£ Integers. 

(10) 

The above Equation is the discrete version of the Hamil- 
tonian for a free particle moving on the lattice. Using 
the Bloch function approach 29] for the periodic system 
in solid state physics, the time independent SE reads 



H\ip„ >= E n \ip n > 



(11) 
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The Bloch states \ip n > can be expanded as a linear com- 
bination of the states \j > localized at node j, 



(12) 



Substituting Eqs. (TO]) and JT2J! into Eq. flTTJ), we obtain 
the eigenvalues (or energy) of the system, 



classical transport where there is a uniform probability 
(1/JV) to find the walker at every node. The time limiting 
probabilities depend on the degeneracies of the eigenval- 
ues, which result in odd, unexpected patterns of limiting 
probability distributions. 



III. TIME EVOLUTION OF THE 
PROBABILITIES 



E n = 2m — 2 cos(Jd n ) 



(13) 



The periodic boundary condition for the network requires 
that the projection of on the state |JV+1 > equals to that 
on the state |1 >, thus 6 n = 2mr/N with n integer and 
n e [0, JV). Replacing \q n > by the Bloch states \ip n > 
in Eqs. I©, 0, we can get the classical and quantum 
transition probability 



n k ,j(t) = \a k ,j(t)\' 



e -tE ne -i(k-j)2^ 



(14) 



N 



j2 z_-m,i 



-it(E„ 



-B,) e -i(fc-i)(n-0 ( 15 ) 



For infinite networks, i.e. 
translates to 



Jim Pk,j(t) 



-Irat 



N^oo 



2tt 



jv 



lim 7T fe j (t) = | 



1 

2^ 



JV -> oo, Eqs. (HI]) and (fJS]) 

e -i9(k-j) e 2t J2?=i cos i d0 

(16) 

-i0(k-j) e 2it T.T=i cos 3 e d9\ 2 



(17) 

Particularly, when m = 1, the network corresponds to a 
cycle graph where each node has exact two nearest neigh- 
bors. The limiting transition probability can be rewrit- 
ten as linijv^oo Pk,j (t) — e" 2mt BesselJ(k — j, 2t) and 
limjv— >oo ^k,3 (t) — \BesselJ(k — 7, 2t)] 2 , where BesselJ 
is the Bessel function of the first kind [3l| . This is consis- 
tent with the result in Ref. [23j. The difference between 
finite and infinite network is that, for infinite networks 
the interference of quantum transport is weak compared 
to finite networks. For larger value of m, the above ana- 
lytical expression could not be further simplified. We can 
calculate the transition probabilities straightly using the 
integrations for the infinite networks. We will show that 
there is some difference of the transition probabilities be- 
tween finite and infinite networks at long time scale. 

Finally, the long time averaged probability between 
two nodes yields 



i(k-j)(n-l) 



(18) 



Interestingly, the long time averaged probability is re- 
lated to the spectral of the networks, in contrast to the 



In this section, we analyze the time dependent prob- 
abilities of the theoretical calculations. The numerical 
determination of the eigenvalues, eigenvectors and inte- 
gration is done using the software Mathematica. Specif- 
ically, we perform our calculations on infinite and finite 
(JV = 100) networks with different connectivity m. 



A. Return probabilities 

The probability to be still or again at the initial node is 
a good measure to quantify the efficiency of the transport 
[32l |. Classically, according to Eq. (fl4| . the probability 
of being at the original node j is 



JV^' 



-tE n 



(19) 



Which only depends on the eigenvalues. The quantum- 
mechanical probability of finding the walker at the initial 
node is given by Eq. (p~5|) . 



(20) 



Which also dependents on the eigenvalues of the system. 
The return probability is independent on the position of 
the initial excitation nodes because of the symmetry of 
network topology. Analogously, employing the relation 
of k = j, we can calculate the return probabilities on 
infinite networks according to Eqs (Til))) and (TTT]) . 

FigQ] shows the return probabilities for CTRWs and 
CTQWs. Consider a CTRW on a network of size N = 
100 and assume the initial excitation starts at node 1. 
FigfT] (a) depicts the temporal behavior of return prob- 
ability with different values of m. There is a power law 
decay (p ~ i~°- 5 ) at the beginning of the transport, but 
after some time p reaches a constant value. This time is 
determined by the time when reaches the equiparti- 
tioned probability 1/JV. The time becomes small when 
the parameter m increases, this indicates that it takes 
less time for the return probability to reach the equipar- 
titioned probability on networks with high connectivity. 
FigfT] (b), (c) and (d) shows the quantum mechanical 
return probabilities for m — 1, m = 2, and m = 3, 
respectively. The dashed curves show the results on net- 
work of JV = 100 and the black solid curves show the 
results on infinite network according to Eq. (|T7|l . The 
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FIG. 1: (Color online) Evolution of the probability of being 
at the initial node 1. (a)Classical return probability pi.i on 
networks of N = 100 with different values of m. pi.i ap- 
proaches the equipartitioned probability 1 /N quickly on net- 
works with high connectivity, (b), (c) and (d) show the evo- 
lution of quantum mechanical return probabilities 7Ti,i with 
m = 1, m = 2 and m = 3, respectively. The dashed curves 
are results on network of size N = 100 according to Eg. (1201) . 
the solid curves are the corresponding results on infinite net- 
works according to Eq. (|17|l . The dashed lines show the scaling 
behavior 7Ti,i ~ t~ . 



FIG. 2: (Color online) Evolution of the probability of find- 
ing the walker at the opposite node 51 when the initial node 
is 1. (a) The classical transition probability P51.1 for infi- 
nite networks and finite network of N = 100 with different 
parameter m. We can see that, the probability of infinite 
network with large connectivity reaches its maximum quicker 
than that of infinite network with small connectivity, (b), (c) 
and (d) are the quantum mechanical transition probabilities 
7i"5i,i for m — 1, m = 2 and m = 3. The solid curves are the 
results on infinite networks, the dashed curves are the results 
on finite networks of N = 100. 



dashed lines indicate the scaling behavior 7ri,i ~ t . 
We note that the return probabilities of finite and infi- 
nite networks agree with each other in small time scales. 
At later times waves propagating on the finite networks 
start to interferer, this leads to the probabilities differ 
and the deviation happens at earlier times on highly con- 
nected networks (with larger value of m). Furthermore, 
the return probabilities oscillate frequently on highly con- 



nected networks and there are more peaks compared to 
networks with small value of m. Such a behavior may be 
attributed to the fact that the interferences on networks 
with high connectivity are stronger than on those with 
small connectivity. 
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m=1 m=2 m=3 



FIG. 3: (Color online) Development of transition probabilities nk,i(t) f° r CTQWs on network of TV = 100 with parameter 
m — 1 (a), m = 2 (b) and m — 3 (c). The initial exciton starts at node 1. 



B. Transition probabilities and transport velocity 

The transition probabilities between two different 
nodes provide us more information about the transport 
process over the whole network. For a finite network 
of N = 100, we consider the probability of finding the 
walker at the opposite node. Fig. [2] shows the transition 
probabilities for CTRWs and CTQWs. In Figure El (a) 
shows the classical transition probabilities psi,i on infi- 
nite and finite networks of size N = 100 with different 
values of m. As we can see, the transition probabilities 
on finite networks with more connectivity approach to 
the equipartitioned probability 1/N quicker than those 
on network with less connectivity. For infinite network, 
the transition probabilities increase with time in the first 
period, and then reach the maxima and decrease in the 
large scale time. Quantum mechanically, the transition 
probabilities for to = 1, m = 2, and m = 3 are shown 
in Fig. [2] (b), (c) and (d). The dashed curves are the 
results for network of N = 100, solid curves are the cor- 
responding results for infinite networks. The transition 
probabilities on infinite networks are smaller than those 
on finite networks at the same time. Interestingly, for 
the same connectivity parameter to, the character time 
t c when the first maximum of the probabilities occur on 
finite networks equals to that on infinite networks, i.e., 
the character time t c is independent on the size of the 
networks. 

The probabilities to go from a starting node to all other 
nodes in time t on a network of size N = 100 with dif- 
ferent values of m are plotted in Fig. [31 The starting 
excitation is located at node 1, and we can see the time 
propagating to the opposite node 51 becomes small on 



networks with large value of to. In addition, the struc- 
ture is quite regular when m = 1. As to increases, the 
pattern becomes irregular. 

In order to compare the transport speed on different 
networks, we define the character time t c as the time 
when the first maximum of the probabilities occur on 
infinite networks. Such definition is held both for the 
classical and quantum transport. For the classical trans- 
port, there is only one maximal value and the character 
time corresponds to the time when the equipartitioned 
probability 1/N is reached on finite networks. Now it 
is natural to ask the question : Does the transport take 
equal time between two nodes of the same shortest path 
length? To address this question, we calculate the tran- 
sition probabilities between two nodes having the same 
value of shortest path length on infinite networks. Fig. [4] 
(a) shows the classical transition probabilities pn.i, P2X,i 
and P3i,i for m = 1, m = 2 and to = 3. The short- 
est path lengths of the two nodes for the three infinite 
networks equal to 10, but the character time t c is small 
for highly connected networks. This indicates that the 
transport is quick on networks with high connectivity 
for CTRWs. For CTQWs, the same conclusion is also 
true, as confirmed by the corresponding plot in Fig.[U(b). 
The character time t c for the quantum transport is much 
smaller than that of the classical one, this supports the 
fact that the quantum walks have efficient quantum al- 
gorithms with respect to their classical counterparts [33[ . 

Fig. [5] shows the character time t c versus the short- 
est path length on networks with different values of to. 
For classical transport (Fig. [3(a)), t c grows faster than 
the shortest path length L. It is found that the rela- 
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FIG. 4: Time evolution of transition probabilities on infi- 
nite networks for CTRWs (a) and CTQWs (b). The initial 
excitation is located at node 1. The solid curves show the 
probabilities of being at node 11 for m — 1 (pn,i in (a) and 
7Tn,i in (b)). Dashed curves show the probabilities of being 
at node 21 for m = 2 (p2i,i in (a) and 7T2i,i in (b)). Dotted 
curves show the probability of being at node 31 for m = 3 
(psi,i hi ( a ) an d 7T31,! in (b)). The shortest path lengths be- 
tween the two nodes are equal, but the time when the first 
maximal value appears are different. 



tionship between the character time t c and the shortest 
path length L can be well described by quadratic equa- 
tion: t c = f3L 2 , where the parameter (3 can be obtained 
by fitting the data. Defining the transport speed v as 
the ratio of L and t c , we find that the classical transport 
speed gets slow for large L while the quantum transport 
speed turns out to be a constant values. We note that the 
transport speed v is large on highly connected networks 
even the two nodes are located at the same distance Li j . 
By fitting the linear relation between t c and L, the quan- 
tum transport velocities are estimated to be about 1.92, 
2.62 and 3.41 for m = 1, m = 2 and m = 3 respec- 
tively. The different behavior of the transport velocities 
between CTRQs and CTQWs is a striking characteristic 




20 j_ 30 



FIG. 5: The chacter time t c as a function of the shortest 
path length L with different values of m for CTRWs (a) and 
CTQWs (b). From the figure, we can see that the classical 
transport gets slower while the quantum transport velocity 
turns out to be invariable for a certain value of m. 



that distinguishes the classical and quantum transport 
processes. 

IV. LONG TIME LIMITING PROBABILITIES 

Now, we consider the long time averaged probabilities. 
Classically, the long time liming probabilities equal to 
the equipartitioned probability 1/iV [23| . Quantum me- 
chanically, the limiting probabilities are determined by 
Eq. (UHJ) but the situation is more complex for different 
network parameters. For m = 1, the spectral (or energy) 
of the system is E n = 2 — 2 cos(6*„), where 9 n — 2mr/N, 
n 6 [0, N). If the network size N is an even number, 
there are two nondegenerate eigenvalues, E N / 2 = 4 and 
Eq = 0, and other eigenvalues have degeneracy 2. The 
limiting probabilities can be written as 



Xij 



2(N - l)/N 2 , 
(N -2)/N 2 , 



if i = j , i - 
Otherwise. 



: j ± N/2, 



(21) 



If the network size N is an odd number, there are one 
nondegenerate eigenvalues En = 0, and the other eigen- 
values have degeneracy 2. The limiting probabilities can 
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FIG. 6: (Color online) The long time averaged probability 
distribution Xk,i f° r CTQWs on networks of size N = 100 
with different values of m. 



be summarized as 



Xij 



_ J {2N-l)/N\ if i=j, 

(N~l)/N 2 , Otherwise. 



(22) 



Which confirms the results in Ref. [30fl. 

For other values of m, the limiting probability distri- 
butions can also be determined according to the degen- 
eracy distribution of the eigenvalues, but such process is 
complicated for large values of m. Here, we report the 
limiting probabilities numerically obtained using the Eq. 
(JTSJ). In Fig. O we display the limiting probabilities on 
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FIG. 7: (Color online) Quantum mechanical limiting proba- 
bilities Xk,i on networks of size N = 75 with different values 
of ra- 



the network of size N = 100 with the starting node 1. 
As we can see, the probabilities for m = 6 and m = 8 
are the same as m = 1 and m = 3. After a careful 
examination, we find that m = 8 and m = 3 have the 
same degeneracy distribution of eigenvalues, and m = 6 
and m — 1 have the same degenerate eigenvalue distribu- 
tion. Particularly, for all the values of m, there is a large 
probability to be still or again at the initial node and at 
the opposite node k — 51. For some values of m, the 
probabilities at the two positions are extremely high, for 
instance, when m — 12, the return probabilities exceed 
0.07. For odd number of network size, there is a higher 
probability to find the walker at the initial node than 
that at other nodes. For networks of size N — 101 and 
m 7^ 50, the limiting probability distribution shows the 
same pattern described as Eq. ([22|) . One may conjecture 
that the pattern of Xh,i does not change when increasing 
parameter m on odd- numbered networks, but this is not 
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FIG. 8: (Color online) (a) Relationship between the quantity 
A(l,50) = (xi,i - X5i,i)/(Xi,i +Xsi,i) an d m on a network 
of size N = 100. The nonzero value of A(l, 50) presents 
asymmetry of probabilities 71-1,1 and 7151,1. (b)A(l, 7V/2) = 

(Xl,i -Xiv/2+i,i)/(xi,i + XJv/2+l,i) versus network size N for 
different values of m. The solid line indicates the power law 
decay A(l,7V/2) ~ N' 1 . 



true for some values of network size N. For instance, on 
networks of size N — 75 with some particular values of 
m, the limiting probability distribution differs from the 
pattern of Eq. (f2"2")l (See Fig. [7]) . It is interesting to note 
that the patterns of Xfc,i are the same for some values of 
to, this feature can be explained by the identical degen- 
eracy distribution of the eigenvalues for different values 
of to. 

As we have shown, if the network size N is an even 
number, there are high probabilities to find the walker at 
the initial node and the opposite node. For some values of 
to, we find that the probability of being at the initial node 
equals to the probability of being at the opposite node. 
However, for some other values of to, this is not true. In 
Ref. [30( , the authors find asymmetry of the probabilities 
for the starting node and its mirror node, their definition 
of mirror node is based on geometry symmetry of the 
network. In this paper, we define the mirror node i of a 
given node i to be its opposite node, i.e., i = i+N/2. We 
find asymmetry of the probabilities of being at the initial 
node and at the opposite node (mirror node) for some 
particular network parameters N and to. Such asymme- 
try is small and not easy to be observed in Fig. For 
a network of size N — 100 and assuming the initial ex- 



citon starts at node 1, we find that asymmetries occur 
at m = 3, 4, 5, 7, 8, 11, 12, 14, 15, 16, 19, 20,.... The asym- 
metrical limiting probabilities are particularly character- 
ized by the difference between xi,i and X5i,i> thus we 
use the quantity A(l,50) = (xi,i - X5i,i)/(xi,i + Xsi.i) 
to detect the asymmetry of the probabilities. In Fig. [5] 
(a), we present A(l,50) as a function of parameter to. 
There are 29 distinct values of to having asymmetrical 
probabilities, which is indicated by the nonzero value of 
A(l,50). 

To reveal a general dependence of the asymme- 
try on the network parameters, we plot the quantity 
A(l,iV/2) = (xi,i-Xjv/2+i,i)/(xi,i+Xjv/2+i,i) as a func- 
tion of the network size N for different values of to, which 
are shown in Fig. [8] (b). For m = 1, the probabilities are 
symmetrical for all the network size N, thus we only show 
the asymmetry for to = 2, to = 3, and to = 4. We find 
that the points break into several clusters, whereas some 
clusters A(l, N/2) decreases with the network size N as 
a power law: A(l, N/2) - N' 1 . 

Except for the asymmetrical probabilities between the 
initial node and the opposite node (mirror node), we also 
find asymmetrical probabilities between other nodes and 
their mirror nodes. In our calculations, we find that such 
asymmetries can be different from the asymmetry of the 
probability of being at the initial node and being at its 
opposite node. For instance, considering a CTQW on a 
network of size N = 100 and assuming the initial exci- 
tation starts at node 1, there are asymmetries between 
Xi+n,i and X5i+n,i i n is an even number) for some values 
of to. The discrete values of to for different asymmetries 
can differ from each other, depending on the precise value 
of N and to. This situation is even more complex and 
requires a further study. 



V. CONCLUSIONS AND DISCUSSIONS 

In summary, we have studied continuous-time quan- 
tum walks on one-dimension ring lattice of N nodes in 
which each node is connected to its 2to nearest neighbors 
(to on either side). Using the Bloch function approach, 
we calculate transition probabilities between two nodes 
of the lattice, and compare the results with the classical 
counterpart. It is found that the transport of CTQW 
is faster than that of the classical continuous-time ran- 
dom walk. We define the transport velocity as the ratio 
of the shortest path length and spreading time between 
two nodes. For network of a given parameter to, the 
transport of CTRWs gets slow with the increase of the 
shortest distance while the transport of CTQWs spreads 
the network constantly. In the long time limit, depending 
on the network parameters N and to, the limiting prob- 
ability distributions of CTQWs show various patterns. 
When the network size N is an even number, the prob- 
ability of being at the original node differs from that of 
being at the opposite node, which also depends on the 
precise value of parameter to. Asymmetrical probabili- 
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ties between other nodes and their mirror nodes also exist 
for some particular network parameters. 

The asymmetry of the limiting probabilities of be- 
ing at a node and being at its mirror node is a novel 
phenomenon, which does not exist in the cycle graph 
with m = 1. However, we are unable to predict which 
particular parameters of N and m are related to such 
asymmetry. Furthermore, in our calculations, we find a 
large value of the limiting return probability for some 
special network topology, for instance, on a complete 
graph in which each pair of nodes is connected, the 
long time averaged return probabilities equal to Xj,j — 
(N 2 - 2N + 2)/N 2 while the other transition probabili- 



ties are Xk,j = 2/7V 2 (k ^ j). This is a striking feature 
of CTQWs which differs from the classical counterpart. 
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